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ABSTRACT 

We study the properties of simulated high-redshift galaxies using cosmological N- 
body/gasdynamical runs from the Overwhelmingly Large Simulations (OWLS) project. The 
runs contrast several feedback implementations of varying effectiveness; from no-feedback, 
to supernova-driven winds to powerful AGN-driven outflows. These different feedback mod- 
els result in large variations in the abundance and structural properties of bright galaxies at 
z = 2. In agreement with earlier work, models with inefficient or no feedback lead to the 
formation of massive compact galaxies collecting a large fraction (upwards of 50%) of all 
available baryons in each halo. Increasing the efficiency of feedback reduces the baryonic 
mass and increases the size of simulated galaxies. A model that includes supernova-driven 
gas outflows aided by the energetic output of AGNs reduces galaxy masses by roughly a fac- 
tor of ^ 10 compared with the no-feedback case. Other models give results that straddle these 
two extremes. Despite the large differences in galaxy formation efficiency, the net specific 
angular momentum of a galaxy is, on average, roughly half that of its surrounding halo, inde- 
pendent of halo mass (in the range probed) and of the feedback scheme. Feedback thus affects 
the baryonic mass of a galaxy much more severely than its spin. Feedback induces strong 
correlations between angular momentum content and galaxy mass that leave their imprint on 
galaxy scaling relations and morphologies. Encouragingly, we find that galaxy disks are com- 
mon in moderate-feedback runs, making up typically ^ 50% of all galaxies at the centers of 
haloes with virial mass exceeding ^ 10^^ A/0. The size, stellar masses, and circular speeds of 
simulated galaxies formed in such runs have properties in between those of large star-forming 
disks and of compact early-type galaxies at z = 2. Once the detailed abundance and structural 
properties of these rare objects are well established it may be possible to use them to gauge 
the overall efficacy of feedback in the formation of high redshift galaxies. 

Key words: galaxies: haloes - galaxies; formation - galaxies; evolution - galaxies; kinematics 
and dynamics. 



1 INTRODUCTION 

The established paradigm for structure formation offers a clear road 
map for galaxy formation. Primordial fluctuations in the dominant 
cold dark matter (CDM) component of the Universe grow via grav- 
itational instability, sweeping baryons into an evolving hierarchy of 
dark matter haloes that grow through mergers of preexisting units 
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as well as t hrough the accretion of material from the intergalac- 
tic medium jWhite & Rees|[l978h . On galaxy mass scales, baryons 
caught in a halo are able to radiate away the gravitational energy 
gained through the collapse, sink to the center of the halo, and as- 
semble i nto the dense aggregati ons of gas and stars that we call 
galaxies telumenthaletal.ll 19851) . 

The structure and morphology of a galaxy results from the 
complex interplay between the time of collapse, the mode of as- 
sembly, the efficiency of cooling, and the rate of transformation of 
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gas into stars (see, e.g. JSteinmetz & Navarr j2002h . Where cooling 
dominates and outpaces sta r formation, baryons co llect into thin, 
rotationally-supported disks (Fall & Efstathiou '1980'). Stars formed 
in these disks inherit these morphological features, but are vulnera- 
ble to swift transformation into dispersion-sup ported spheroids dur- 
ing subsequent merger events l lToomre|[l977l) . Disks may re-form 
if mergers or accretion bring fresh supplies of cooled gas, making 
morphology a constantly evolving rather than an abiding feature of 
a galaxy ( Cole etal.. 2000 ; Robertson et al. 2006). 

The galaxy formation scenario driven by gravitational col- 
lapse and radiative losses outlined above is compelling, but incom- 
plete. Indeed, cooling is so effective at early times that, unless im- 
peded somehow, most baryons would be turned into stars in early- 
collapsing protogalaxies, which would then merge away to form 
by the present time a majority of spheroid- dominated remnants , 
in ve hement disagreement with observations I White&Ree3l 19781 : 
ICole 1991 ; White & Frenk 1991). The problem is compounded by 
the fact that, during mergers, cooled gas tends to transfer its an- 
gular momentum to the surrounding dark matter halo. As a result, 
even in cases where disks could form, their str uctural properties 
would be at od ds with those of spiral galaxies ( N avarro & Bend 
ll99lUNavarro et al. 1995; Navarro & Steinmetz 199^^ 

A gas heating mechanism that prevents runaway cooling and 
that regulates the formation of stars in step with mergers and ac- 
cretion is widely believed to be the most likely solution to these 
problems. The energetic output from evolving stars and supernovae 
is a natural candidate. It scales directly with star formation and, in 
a typical galaxy, the total energy released by supernovae can be 
comparable to the binding energy of the baryons. Thus, if chan- 
neled properly, feedback energy from supernovae may temper the 
gravitational deposition of cooled gas into a gal axy and effectively 
self-regulate its star formation history dWhite & Frenk. 1991) . 

The even standing of gravity, feedback and cooling may thus 
help reconcile the observed galaxy population with hierarchical 
clustering models, but it comes at the price of complexity: the 
main structural properties of a galaxy, such as stellar mass, ro- 
tation speed, and morphology, are then expected to depend on 
details of its assembly history, such as the exact tim ing, geome- 
try and mass spectrum of accretion events (see, e.g., Abadi et al. 
20033; b; Me za et al. 20()3; Governa to et al. 2007; Zavala et al. 
20081 ; IScannapieco et al.ll2009l : lGovernato et allboiol) . 

Such sensitivity to feedback has held back progress in di- 
rect simulation of the process of galaxy formation. As recent work 
demonstrates, different but plausible implementations of feed- 
back within the same dark halo lead to galaxies of very different 
mass, morphology, d ynamics, and star formation history (see, e.g., 
lOkamoto et afluOOsh . Numerical parameters may thus be tuned to 
reproduce some properties of individual galaxies, but at the expense 
of wider predictability in the modeling. 

These results suggest that further progress in the subject re- 
quires the testing of different feedback schemes on a statistically 
significant sample of dark haloes formed with representative as- 
sembly histories. The viability of each feedback implementation 
may then be assessed by contrasting the statistics of such samples 
with observational constraints such as the stellar mass function, 
clustering, color distribution, and scaling laws. 

We take a step in this direction here by analyzing a sub- 
set of cosmological N-body/gasdynamical simulati ons from the 
Overw helmingly Large Simulations (OWLS) project jSchave et al.l 
l2010h . We present results regarding the morphology, stellar mass, 
and angular momentum content of galaxies assembled at z = 2, 
and compare them with the few observational constraints available 



at that epoch. We limit our analysis to the z — 2 galaxy population 
because most high-resolution OWLS runs follow volumes too small 
to be evolved until z = Q.ln future papers, we plan to extend this 
analysis to the present-day galaxy population using samples drawn 
from the closely-related GIMIC project, designed to follow a few 
representative vol umes selected from the Millennium Simulation 
dCrain et al.ll2009h . 

The paper is organized as follows. In Sec.[2]we present a short 
overview of the simulations and feedback models. We then present 
our main numerical results in Sec.[3]and analyze them in the context 
of available observational constraints in Sec.|4] We end with a brief 
summary in Sec.|5] 



2 THE NUMERICAL SIMULATIONS 
2.1 The OWLS runs 

The OWLS project consists of a suite of ~ 50 different cosmo- 
logical N-body/SPH simulations that follow the evolution of dark 
matter and baryons in boxes of 25 and 100 h^^ Mpc (comoving). 
Each box is run many times, varying the numerical implementation 
of various as pects of the gas coo ling, star formation and feedback 
modules (see ISchave et al.ll20IOl for further details). 

We have selected for our analysis nine 25 /i^^ Mpc-box OWLS 
runs, eight of which explore different feedback implementations 
with 512'' dark matter and 512^ baryonic particles whilst keeping 
other subgrid parameters constant, such as the stellar initial mass 
function (IMF), the star formation threshold and its efficiency. The 
ninth repeats one of the runs, at 8x lower mass resolution (and 
2x lower spatial resolution), in order to provide some guidance 
regarding the sensitivity of our results to numerical resolution. 

All simulations assume a standard WMAP-3 ACDM cos- 
mogony, start at Zi = 127 and, because of their small box size, 
they have only been carried out to z = 2. We adopt this cosmology 
for all physical quantities listed here. We make explicit the depen- 
dence on the Hubble constant, h, for simulation parameters, but 
drop the h dependence and adopt h = 0.73 when comparing with 
observations. 

The high-resolution runs have a comoving gravitational soft- 
ening lengthscale of 1/25 of the initial mean inter-particle spac- 
ing at high redshift. These are switched later to a fixed physical 
value so that the softening never exceeds 0.5 kpc (physical). 
The mass per baryonic particle is ~ 1.4 x l(f h~^yiQ and 4.5 
times higher for the dark matter component. All runs assume that 
the Universe is reionized at z = 9 (for H) and at z = 3.5 (He) by a 
bath of energetic photon s whose properties evolve as proposed by 
iHaardt&MadaulhoOlh . 

Table[T]summarizes the most important numerical parameters 
of the simulations, as well as the cosmological parameters. 



2.2 Subgrid gas physics 

Baryons are assumed to trace the dark matter distribution at the 
initial redshift. Whilst in gaseous form, they are followed hydrody- 
namically and are subject to pressure gradients and shocks. Radia- 
tive cooling and heating is implemented following Wiersma et al. 
(2009a), which also accounts for the photo-ionisation of metals due 
to the UV background. 

In collapsed structures, gas can cool and sink to the cen- 
ter of these haloes, where it may reach high overdensities before 
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Table 1. Simulation parameters 



fiCDM 

Ha 

h 
n 

Reionization redshift 
Mass per DM particle 
Mass per baryonic particle 
Number of particles 
Box size 



0.238 
0.1962 
0.0418 
0.762 
0.74 
0.73 
0.951 
9 (H), 3.5 (He) 
= 6.3 X lO^/i-iMg 
= 1.4 X lO^h-iM 
2 X 512^ 
25 h-^ Mpc 







2.4.1 Thermal Feedback 

The simplest possibility, which we label "thermal feedback" (ThF), 
is to use the supernova energy to raise the internal e nergy of th e 
surrounding gas particles. As reported in earlier work jKatzlll992h . 
these regions typically have such short cooling times that the in- 
jected energy is quickly radiated away, with little hydrodynami- 
cal effect on the surrounding gas. As a result, thermal feedback 
is rather inefficient, and has little effect in regulating gas cooling 
and star formation, even though the implementation here follows 
the stochastic heating method described in S chaye et a l. ( 2010) and 
presented in more detail in Dalla Vecchia et al. (in preparation), 
which is more resilient to numerical resolution li mitations than th e 
implementations adopted in earlier work (see also lKav et al .l2003h . 



turning into stars. With limited numbers of particles, these re- 
gions are poorly resolved and vulnerable to numerical instabili- 
ties, such as artificial clumpi ng and fragmentation. As discussed by 
ISpringel & Hemguistl ( |2003|) . these shortcomings can be alleviated 
by adopting, in high-density regions, a multi-phase description for 
the gas where the effective equation of state differs from the simple 
ideal gas law. In practice, we impose a polytropic equation of state 
(PEOS; P oc p^, with 7 = 4/3) on all gas particles whose den- 
sity exceeds a critical value of Uc — 0.1 cm""', the density above 
which the gas is expect ed to be multiphase and unstable to star 
formation jSchavell2004h . This choice ensures that the Jeans mass 
in high-density regions is independent of p, effectively suppress- 
ing artificial clumping and reducing the d ependence of star forma- 
tion a lgorithms on numerical resolution l iSchave & Dalla Vecchial 
l2008i) . 

2.3 Star formation algorithm 

Star formation is implemen ted as described in detail by 
IScha ve & Dall a Vecchij j2008[) . In brief, stars form out of PEOS 
gas particles with pressure-dependent parameter s chosen to repro - 
duce a Kennicutt-Schmidt law with index 1.4 (lKennicuttlll998h . 
We assume a Chabrier initial mass function jChabrieill200'3l) m or- 
der to take into account the enrichment and energy injected into the 
surroundings of young star particles by the explosion of SNII and 
SNIa supemovae. The energy per supernova explosion is chosen 
to be 10^^ ergs. These events, together with mass loss from inter- 
mediate mass stars, pollute neighboring gas particles with metals, 
as described in Wiersma et al. (2009b). We track 1 1 species and 
include them in the computation of the cooling function follow- 
ing (Wiersma et al 2009a) in all our runs, with the exception of 
the "NoF" model described below. For the latter case, chemical en- 
richment is modelled in the same way but is not considered in the 
computation of cooling, which instead assumes primordial abun- 
dances. 

2.4 Feedback Models 

The runs we analyze here explore alternative feedback implementa- 
tions where the total amount of energy injected by supernovae into 
the surrounding ISM is kept constant, but the numerical algorithm 
used to inject this energy is varied. All runs that include feedback 
from core collapse supernova feedback assume a total energy input 
of 10^^ ergs per solar mass of stars formed, 40% of which is in- 
vested into driving outflowing winds. The remainder is assumed to 
be lost to radiative processes. 



2.4.2 Kinetic Feedback 

A second possibility is to invest part of the feedback energy directly 
into gas bulk motions, with the aim of allowing gas to outflow from 
regions of active star formation, thus increasing the overall effi- 
ciency of feedback. These wind models are characterized by a cou- 
ple of parameters: a "mass loading" factor, 77, specifying the num- 
ber of gas particles amongst which the injected energy is shared, 
and a "wind velocity", Uw, characterizing the kinetic energy of the 
outflow. For given energy, -q and are related by a constant -q . 

For the reference model in OWLS (WF2 in our notation) the 
wind velocity = 600 km/s is c hosen partly motivate d by obser- 
vation of local starburst galaxies jVeilleux et J]|2005}) . The mass 
loading t] is thus fixed to 2 particles as the integer number that best 
reproduces the peak in the cosmic star formation history. This com- 
bination of Uw and r) imply that 40% of the total energy liberated 
by supemovae impacts the kinematics of the surrounding gaseous 
medium. Because all these parameters are highly uncertain, in the 
OWLS runs we contrast the results obtained with three different 
values of 77=1,2,4 (we refer to these runs as WFl, WF2, and WF4, 
respectively). The wind velocities in each model are adjusted so 
that the same amount of energy (40%) is input in every case (see 
Table|2l. Further details can be found in Dalla Vecchia &. Schaye 
(2008) and Schaye et al. (2010). WF2LR is equivalent to WF2 but 
run at 8x poorer ma ss resolution and 2 x poorer s patial resolution. 

As discussed bv lSpringel & Hemguistl j2003l) . a possible mod- 
ification that can enhance feedback efficiency is to temporarily "de- 
couple" the wind particle(s) hydrodynamically from the surround- 
ing ISM. This facilitates large-scale galactic outflows and regulates 
star formation more effectively by enha ncing the removal of gas 
from active star-forming regions (see e.g..l Dalla Vecchia & Schava 
[2008). The criterion for re-coupling particles to the gas is as in 
Springel & Hernquist, and occurs as soon as either: a) the density 
has fallen to 0.1 Wc, where ric is the density threshold for star for- 
mation, or 6) 50 Myr have elapsed since decoupling. We label this 
run WF2Dec. 

A further run probes the possibility that the efficiency of feed- 
back should correlate with the local density of the gas. We explore 
a model in which the wind velocity and mass loading are related to 
the gas density by Hw oc p^''*' and rj oc p"^^'^ . This guarantees that 
the wind velocity scales with the local gas sound speed {v„ oc Cs) 
given the aforementioned effective PEOS that holds in star-forming 
regions. The and rj relations are normalized so that, at the gas 
density corresponding to the star formation threshold (ric = 0.1 
cm^'^), they match — 600 km/s and rj = 2 particles, consistent 
with the WF2 run. We will refer to this run as WDENS. 
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2.4.3 AGN Feedback 

Our next model enhances feedback by adding to the WF2 feed- 
back the extra energetic input from AGN. This model, which we 
refer to a s AGN, for short, follows the numerical procedure intro- 
duced bv lBooth & S chave (2009) and summarized in Schave et al. 
1 I2OIOI) . Seed black holes with mass msecd = 9 x 10*i\fo are 
placed at the center of all haloes that exceed a threshold virial 
mass of 4 X lQ^'^h~^MQ. BHs can then grow by mass accre- 
tion and mergers with other BHs. Gas accretion onto the BH is 
modelled according to a modified version of the Bon di-Hoyle- 
Lyttleton l lBondi & Hovlell944lHovle & Lvttletonll 19391) formula: 
w.accr = a 47r JMbh pI (Cs + v'^)^^'^ , where rriBH is the mass 
of the BH, V is the velocity of the BH with respect to the ambient 
medium, Cs is the local speed of sound and p the local gas den- 
sity, a is an extra "efficiency" parameter that did not appear in the 
ori ginal versions of the Bondi-Hoyle formula but was introduced 
by Springel et al.l ( |2005|) . who set it to alpha = 100, to account for 
the finite numerical resolution and for the fact that the cold, inter- 
stellar phase is not explicitly modeled. In our model this factor is 
set to unity in the regime that the physics is modeled correctly, but 
increases with the ISM density (i.e., for particles on the PEOS, see 
Booth & Schaye (2009) for further details and discussion). 

A fraction e/ of the total radiated energy due to the mass ac- 
cretion onto the BHs is assumed to couple to the surrounding ISM. 
This efficiency is set to = 0.15 to match local constraints on 
the n umber density (see, e.g., Ishankar et alj|20()4l : iMarconi et al.l 
I2OO4') as well as relations between BHs and host galaxy properties 
(Haring & Rix 2004; Tremaine et al. 2002), both at redshift zero. 
AGN feedback is implemented as a thermal injection of energy (as 
opposed to the kinetic prescri ption used to model the stellar feed- 
back), in the way described in iBooth & Sch"av3 ( |2009|) . Because it 
combines the supernova and AGN energetic outputs, the AGN run 
is the most effective feedback model tried in our series. 



2.4.4 No Feedback 

Finally, mainly for comparison purposes we also analyze a run that 
follows star formation like in the other implementations but ne- 
glects all energy injection into the ISM due to either supemovae 
or AGN. Gas cooling in this "no feedback" model, NoF, adopts the 
cooling function of a gas with primordial abundances, but in the ab- 
sence of feedback this is only a minor difference that has little im- 
pact on the results. The NoF model stands at the opposite extreme 
as AGN, allowing for unimpeded transformation of gas into stars in 
regions able to collapse and condense into galaxies. Although un- 
realistic as a galaxy formation model, it serves to provide a useful 
framework where the relative importance of feedback effects may 
be gauged and understood. 

Table [2] summarizes the relevant parameters of each feedback 
implementation. For ease of reference, we also quote in each case, 
the name used to label each simulation bv lSchave et alj ( I2OIOI) . 



Table 2. Parameters of the different feedback modefs probed in each run. 
First and second columns list t he short name (us ed throughout this paper) 
and the name originally used in lSchave et aljliol o^. respectively. The third 
and fourth columns list the mass loading (77) and wind velocity (t)w) pa- 
rameters of each model. The WF2Dec is the only model where wind gas 
particles are temporarily kinematically decoupled from the surrounding gas. 
This aids the removal of gas from galaxies and results in increased feedback 
efficiency. The characteristic density ?iw used for scaling the WDENS wind 
parameters is that corresponding to the star formation threshold: ric = 0.1 
cm~3. 



Short name OWLS name r] t)w 

[particles] [km/s] 



NoF NOSNJSfOZCOOL 

ThF WTHERMAL 

WF4 WML4V424 4 424 

WF2 REF 2 600 

WFl WML1V848 1 848 

WF2Dec WHYDRODEC 2 600 

WDENS WDENS 2(n/n„)-i/^ 600(n/nw)^/® 

AGN AGN 2 600 



3 NUMERICAL RESULTS 

3.1 The halo sample 

Our sample consists of all galaxies at the centers of haloes with 
virial mas £]A/vir > 10" /i"' Mq. There are about 150 haloes at 
z = 2 in each 25 h^^ Mpc-box OWLS run with masses between 
10" Mq < IVUir < 3 X lO^^ft"^ Mq. The median of the 
sample is Mvir ~ 1-8 x 10^^ Mq. Ha los are identified by 
the substructure fi nding algorithm SUBFIND dSpringel et allbOOll: 
iDolag et al.ll2009l) and centers are defined by the minimum of the 
potential. All runs use the same initial conditions, and therefore the 
number (and identity) of haloes selected for analysis is roughly the 
same in each simulation. 

We begin with an overview of the properties of the gaseous 
component within the virial radius and its halo mass dependence 
(Sec. l3.2t . and follow on with a description of the properties of the 
stellar component of the central galaxy (Sec. |3.3l l. We discuss the 
link between feedback and morphology in Sec. 13.41 and compare 
the number of massive galaxies in various runs in Sec. 13. 51 We end 
this section by discussing the mass and angular momentum of cen- 
tral galaxies, as well as their dependence on feedback (Secs. [T6l and 
13. 7t . before proceeding to compare these results with observations. 

3.2 Gas within the virial radius 

Fig. [T] illustrates the distribution of gas within the virial radius in 
four haloes selected from the WF2 run al z — 2. Each panel 
corresponds to haloes differing by consecutive factors of two in 
virial mass. The box size in each panel has been adjusted to the 
virial radius of each halo. Only gas particles within the virial ra- 
dius are shown, and have been colored according to their den- 
sity/temperature. 



^ Virial values are measured at or within the virial radius, Tvir, of a halo, 
defined as the radius where the mean inner density exceeds the critical den- 
sity of the universe by a factor Avir(z) = IStt^ + 82f{z) - 39/(2)^. 
Here f{z) = [ ^Mfl + z)^/(nM(l + zf + Qa))] - 1 and = 
f^CDM + f^bar <Brvan & Normanll 19981) . A^ir ~ 170 at 2; = 2 for our 
choice of cosmological parameters. 
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Figure 1. Gas particles within the virial radius of four WF2 haloes spanning the mass range of systems selected for analysis. The virial mass and radius ai'e 
given in the label of each panel. Gas particles are colored according to temperature: red, green, and blue correspond to particles in the hot, warm, and cold 
phases, respectively. Hot particles are those with T > (l/4)Tvir, where T^ir = 35.9(yvir/ km s~^)^ K is the virial temperature of the halo. Cold particles 
are those with T < 3 X 10^ K. Warm are those with intermediate temperatures. Cyan particles denote dense, star-forming gas in the PEOS phase. Particles 
are plotted sequentially in order of descending temperature, so colder particles may occult hotter ones in regions of high density. Small circles show the radius, 
'■gal = 0-15 ''vir, used to define the central galaxy. 



Red particles are those with temperatures exceeding 
(1/4) Tvir, where Tvir = 35.9 (Kir /km s'^)^ K is the virial tem- 
perature of a halo (Vvir is the circular velocity at rvir). Gas particles 
in this "hot phase" are all found in a low-density, largely pressure- 
supported atmosphere that fills the halo out to the virial radius. The 
virial temperature is ~ 10^ K for haloes with Vvir ~ 170 km/s, 
about the median virial velocity range spanned by our sample. 

The fraction of gas in the hot phase increases with halo mass; 
it makes up 68% of all the gas within rvir in the most massive halo 
but only 21% in the least massive system shown in Fig.[i|. This is 
a result of the steady increase in cooling time with increasing halo 
mass, which favors the formation of a hot tenuous gas atmosphere 
in massive systems. 

Particles in green are those in the "warm" phase, which we 
define as those satisfying 3 x 10*K< T < (1/4) Tvir. These are 
particles at moderate overdensities, and make up a small fraction 
of all the gas within rvir; from ~ 7% in the most massive halo to 
^ 15% in the least massive one. This gas typically traces mate- 
rial accreted relatively recently, which has yet to be pressurized by 



shocks, or material ejected during accretion events in "tidal tails" 
that expand and cool as they recede from the center. Because ac- 
cretion occurs frequently through filaments, and tidal tails are like- 
wise highly asymmetric, the warm component distribution is non- 
uniform, with discernible large-scale features suggestive of recent 
mergers and accretion events. 

Cold (T < 3 X lO" K) gas of moderate density (n < ric = 0.1 
cm"'') is shown in blue, and is rather clumpy in appearance. Large- 
scale features similar to those noted for the warm component are 
also visible here, suggesting that this is also mostly gas recently ac- 
creted or affected by accretion events. In terms of mass, this com- 
ponent is negligible (~ 5%) in the 3.7 x 10^^ h^^ Mq halo but in- 
creases in importance with decreasing halo mass. Indeed, it makes 
up ~ 30% of all the gas in the 3.1 x lO" Mq halo shown in Fig.^ 

The star-forming gaseous component is, by definition, the 
densest (n > Uc = 0.1 cm~'^), and is shown in cyan in Fig. [T] 
Most of this gas is at the bottom of the potential well of the main 
halo and of its substructure haloes, and makes up between 20 and 
30% of the gas within rvir, with little dependence on halo mass. 
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Figure 2. Zoomed-in view of the galaxies at the centers of the haloes shown in Fig.[T] Colors are as described in the caption of that figure, except that yellow 
now denotes "star" particles. The circles show the galaxy "radius", Tg^i = 0.15 rvii-. Each box has been rotated so that the spin axis of the PEOS gas is 
aligned with the z axis of each panel. This "edge-on" projection emphasizes the presence of disk-like structures in all four haloes. Besides r'gai, labels in each 
panel specify the baryonic mass of the galaxy, Mg^i, and the spin parameter of the surrounding halo, A. 



The generally strong halo mass dependence of the various 
gaseous phases highlights the different modes of accretion that 
shape the evolution of a central galaxy. In massive haloes galaxies 
grow by accreting cooled material from the surrounding reservoir 
of hot gas, whereas in low mass haloes the gas is likely to flow 
virtually unimpeded to the central regions, where it ma y be swiftly 
accre t ed into the central galaxy (White & Frenk 199jJ; 



2005 
200S 



Dekel & Birnboin l200d .Bimboim et aL,2007l : 



Keres et al 
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Brooks et al.ll200S ). These different accretion modes high 



light the complex assembly history of a galaxy, a complexity that 
is further compounded by the effects of feedback that we discuss 
below. 



halo. As is clear from the figure, this definition includes virtually 
all stars and dense gas obviously associated with the galaxy. 



It also emphasizes the halo mass dependence of the various 
phases in which baryons may flow into the central galaxy. As dis- 
cussed above, whereas galaxies in low mass haloes grow through 
the smooth accretion of cold gas, a fair fraction of the star forming 
gas in the most massive systems include "clouds" that condense out 
of the hot and warm phases. Little star formation happens in these 
clouds, however, since their typical densities are well below those 
reached in the main body of the galaxy. 



3.3 Central galaxies 

Fig.|2]shows a zoomed-in view of the four WF2 haloes depicted in 
Fig-El including the stellar component, which is shown in yellow. 
The circle centered on the main galaxy indicates the radius, rgai = 
0.15r-vir, that we use to define the central galaxy inhabiting each 



Gas turns swiftly into stars once it settles into a dense, thin, 
rotationally supported disk in the central galaxy. In systems that 
avoid major mergers, the stellar component inherits the disk-like 
structure of the gaseous component. All 4 galaxies shown in Fig.|2] 
sport well-defined stellar disks, which have been rotated to be seen 
"edge-on" in this figure. Disks of gas and stars are indeed quite 
common in the WF2 run that we have chosen to illustrate the main 
general features of our simulated galaxies. 
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Figure 3. Face-on view of the central galaxy formed at the center of an Mvir = 1.2 X 10^^ h^^ Mq halo. All panels con'espond to the same halo, but in runs 
with different feedback implementations, as labelled in the bottom right of each panel. Only baryonic particles are shown. Colors indicate gas temperature, 
classified as hot (red), warm (green), cold (blue) and star-forming (cyan). See the caption to Fig.[T]for details. Yellow dots correspond to "star" particles. 
The circle in each panel indicate the radius used to define the central galaxy, rgai. Each galaxy has been rotated so that it is seen "face on"; i.e., the angular 
momentum of the PEOS gas is aligned with the line of sight of the projection. The mass in stars and gas within rg^i is labelled in each panel (units are 






Figure 4. Same as Fig. [5] but each galaxy has been rotated so that it is seen 
component in ordered rotation. 



'edge-on". Labels in each panel give the fraction of kinetic energy of the stellar 



3.4 Feedback and morphology 

Varying the feedback implementation has a dramatic effect on the 
properties of central galaxies. We illustrate this in Figs. [5] and |4l 
where we show, for the same dark matter halo, how the appearance 
of its central galaxy varies with feedback. Although the assembly 



history of the dark halo is identical in all cases, differences in feed- 
back lead to drastic variations in the stellar mass, gaseous content, 
and morphology of the central galaxy. 

When feedback is inefficient, such as in the ThF and WF4 
runs, a stellar disk is clearly present, but its mass is small com- 
pared with that of the spheroidal component. This is because most 
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NoF ThF WF4 WF2 WFl WF2Dec WDENS AGN 



Figure 5. Number density of galaxies with stellar mass exceeding 5 X 
10^ Mq, shown for the vaiious feedback implementations explored in this 
paper. Runs ai'e labeled as in Table|2] and ai'e listed in the abscissa roughly 
in order of increasing feedback efficiency. The open circle correspond to 
the WF2 low resolution run WF2LR. "Error bars" denote sfN uncertain- 
ties corresponding to the number of systems in our computational box. The 
shaded area outline observational constraints fro m estimates of the galaxy 
stellar mass function at z = 2, as compiled by iMarchesini et al] j2009l) . 
Note the strong decline in the number of massive galaxies as a function of 
increasing feedback efficiency. 



stars form in early collapsing protogalaxies which are later stirred 
into a spheroidal component when these subsystems coalesce to 
form the final galaxy. The extreme case is NoF, where the absence 
of feedback allows for early and highly efficient star formation that 
converts most of the available gas into stars. The large number of 
satellites seen around the NoF central galaxy is also a result of the 
lack of feedback. This preserves star formation even in small sub- 
haloes, where modestly energetic feedback might lead to drastic 
changes in the availability of star formation fuel and in the total 
mass of stars formed. 

When feedback effects are strong, such as in the WF2Dec, 
WDENS, and AGN runs, fewer stars form since the gas is con- 
stantly pushed out of star-forming galaxies by outflows. These out- 
flows also disrupt the smooth settling of gas into disks and its grad- 
ual transformation into stars. In the most extreme case (AGN), the 
gas outflows are so violent that there is little gas left in the central 
galaxy. In none of these cases do central galaxies have an extended 
and easily recognizable stellar disk component. 

As may be seen from Fig. [4] more moderate feedback im- 
plementations, such as WF2 and WFl, yield systems with a well- 
defined stellar disk, and a gas/stellar mass fraction of roughly 1:1. 

This impression is corroborated quantitatively by the fraction 
of stellar kinetic energy in ordered rotation: 



K,^t/K- with K,ot = ^(l/2)m(j./R)'^ 



(1) 



Here, m is the mass of a star particle; is the z-component of 
the specific angular momentum, assuming that the z-axis coincides 
with the angular momentum vector of the galaxy, and R is the 



(cylindrical) distance to the z-axis. n^ot is listed in each panel of 
Fig- El for the stellar component: it is highest for WF2, and mini- 
mum for AGN. 



3.5 Feedback and massive galaxies 

A robust way of assessing the effectiveness of the various feedback 
implementations explored in these runs is to compute the abun- 
dance of massive galaxies that each predicts. Because the total 
amount of stars formed decreases as the feedback becomes more 
effective, the abundance of massive galaxies is expected to depend 
sensitively on feedback. This is shown in Fig. [5] where we plot, 
for each implementation, the number of galaxies (per unit volume) 
with stellar masses exceeding 5 x 10^ Mq . This mass threshold is 
chosen to roughly coincide with 10, 000 baryonic particles 

The runs in Fig. [5] are ranked, from left to right, in order of 
decreasing number of massive galaxies (i.e., increasing feedback 
efficiency). This figure confirms the strong effect of feedback on 
the abundance of massive galaxies. For example, the AGN run has 
~ 16 times fewer such galaxies than the run with only thermal 
feedback, ThF, and ~ 60 fewer than NoF, the model without feed- 
back energy injection. 

Not only does the total amount of feedback energy matter, but 
also the manner in which it is injected. Indeed, large differences are 
also obtained for models with the same feedback strength (as mea- 
sured by the total feedback energy per unit stellar mass formed), 
but that differ in the combinations of mass loading and wind ve- 
locities: ThF, WF4, WF2, WFl, WF2Dec and WDENS all assume 
that 40% of the available supernova energy is invested into winds, 
yet their predictions for the number of bright galaxies differ by up 
to a factor of ~ 4. 

The results do not seem to depend dramatically on numeri- 
cal resolution, as shown by the good agreement between the WF2 
and WF2LR runs (the latter is shown with an open symbol in 
Fig.[5]l- Reducing the number of particles by a factor of eight (as in 
WF2LR) brings down the number of A/gai > 5 x 10^ Mq galaxies 
in the box from 133 to 123. The trends shown in Fig.[5]are therefore 
unlikely to be an artifact of limited numerical resolution. 

The shaded band in Fig. [5] indicates the expecte d number of 
massive galaxies, taken from IMarchesini et al.l ( |2009|) . after inter- 
polating their fits to 2 = 2 and correcting volume elements to ac- 
count for the different cosmology assumed in their work. The band 
aims to represent uncertainties due to photometric redshift inac- 
curacies, cosmic variance and systematics in the modeling. Notice 
that the bright end of the luminosity function traces the abundance 
of the most massive objects present at 2; = 2 and is, as such, par- 
ticularly sensitive to the adopted cosmological parameters. 

Given these large uncertainties, it would be premature to use 
Fig.|5]to rule in or out any particular implementation of feedback 
but, as the data improve, it might be useful to revisit this issue to 
learn which feedback modeling procedure is favored or disfavored 
by the data. For clarity, many of the plots in the analysis that follows 
will focus on 4 cases that span the full range of feedback strength 
shown in Fig.lH i.e., NoE WF2, WF2Dec, and AGN. 

We end by noting that different observational diagnostics, 
such as the specific star formation rate or the gas content as a func- 
tion of galaxy luminosity, could be used to provide further con- 
straints on the viability of each feedback model. We plan to present 
a detailed analysis along these lines in a future paper (see Haas et 
al., in preparation). 



Feedback and the structure of galaxies at z=2 9 




Figure 6. Left: Central galaxy mass as a function of viiial mass for all haloes selected in the WF2 ran. The top dashed line indicates the mass in baryons in 
each halo corresponding to the universal baryon fraction, /bar = ^b/^M = 0.175, adopted in the simulations. Dots correspond to the baryon mass of the 
central galaxy (i.e.; within rg^i); the thick solid (blue) curve tracks its median as a function of A/vir- The two bottom curves track the median for the stellar 
(red, dot-dashed) and gaseous (green, dashed) mass within Tgai. The thick dotted (magenta) line shows the median of the total mass in baryons within the 
virial radius, r^i^. Right: Galaxy formation "efficiency", rjgai = Afgai/(/bar^^vir)> as a function of halo virial mass for the various runs. For clarity, only 
the results coiTesponding to 4 selected rans are shown, spanning the effective range of feedback strength, from the "no feedback" (Nop) case to the AGN case, 
where feedback effects are maximal. Cases not shown fall between these two extremes. The scatter around each curve is large, typically ^ 0.19 dex rms. 
The thin line labelled WF2LR has the same physics are WF2 but 8 X poorer mass resolution and 2 X lower spatial resolution. Note that the galaxy formation 
efficiency is, on average, very sensitive to feedback, but only weakly dependent on halo mass, at least for the range of masses considered here. 



3.6 Galaxy masses 

The stellar and gaseous masses of galaxies assembled at the centers 
of dark matter haloes are determined largely by the virial masses of 
the systems, modulated by the efficiency of radiative cooling and 
the regulating effects of feedback. We show this in the left panel 
of Fig. [6] for all galaxies selected from the WF2 run. The dots in 
the figure correspond to Afgai, the total baryon mass within the 
radius, rgai, used to define the central galaxy; the solid line traces 
the median as a function of Afvir. As expected, the central galaxy 
mass correlates well with Afvir, albeit with fairly large scatter (the 
global rms about the median trend is ^ 0.19 dex). 

The top dashed line in this panel indicates the mass, /barMvir, 
galaxies would have if all baryons in the halo have assembled at the 
center (the universal baryon fraction is /bar ~ Qb/i^M ~ 0.175). 
The thick dotted magenta line shows the median baryon mass 
within rvir as a function of halo mass. This shows that massive sys- 
tems have retained all baryons within the virial radius, but also that 
the effects of feedback are clear at the low mass end: lO" h'^ Mq 
haloes have only retained about half of their baryons within the 
virial radius. Of those, only one third or so have collected in the 
central galaxy. 

Thus, the "efficiency" of galaxy formation, as measured by the 
mass of the galaxy expressed in units of the total baryon mass corre- 
sponding to its halo, r^gai = Mgai/(/bar A/vir), iucrcases steadily 
with halo mass, from ~ 10% in 10^^ h~^MQ haloes to a maxi- 
mum of roughly 40% for A/vir ~ 5 x lO" H'^Mq. There is also 
indication that the efficiency decreases in more massive systems, to 
roughly ~ 30% in the most massive haloes. 

These trends (i.e., low galaxy formation efficiency in low 



and high mass haloes) are qualitatively in line with what is re- 
quired to reconcile the shape of the galaxy lu minosity function 
with the dark matter halo mass function (see , e.g., lYang et alj20o3 ; 
IConrov & Wechsleij|2009l ; lGuo et al.ll2010l) . Feedback is the main 
mechanism responsible for reducing efficiency in low-mass haloes. 
Together with long cooling timescales, it also helps prevent the for- 
mation of too massive galaxies in high-mass haloes. 

Although the trends seem qualitatively correct, it remains to be 
seen whether a model like WF2, evolved to 2: = 0, is able to satisfy 
the stringent constraints placed by the stellar mass function in the 
local Universe. Indeed, the recent estimate of iLi & White! i2009h 
suggest that only 3.5 % of all baryon s in th e Universe are today 
locked up in stars, and McCarthv et al. argue that supernova 

feedback alone is not enough to ensure such a low efficiency of 
transformation of baryons into stars. Since the runs we analyze here 
have only been evolved to z = 2, we are unable to address this issue 
in a conclusive manner, but we plan to return to it when extending 
the present analysis to the GIMIC simulations. 



3. 6. 1 Feedback dependence 

The right panel of Fig.|6]shows that the overall efficiency of galaxy 
formation is quite sensitive to feedback. Each curve here tracks the 
median trend of r^gai with A/vir for different runs. As expected from 
the discussion of Fig. [5] r/gai is highest for NoF and lowest for 
AGN, with more moderate results for WF2 and WF2Dec, as well 
as the other runs, which are omitted from this panel for clarity. 

Clearly, not only the total feedback energy input, but also the 
details of its implementation can affect dramatically the galaxy for- 
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mation efficiency. Central galaxies in the NoF case can be up to 
10 X more massive than in the AGN run. WF2 galaxies are a fac- 
tor of two to three more massive than those formed in WF2Dec, 
although the only difference between these two runs is the choice 
to "decouple" hydrodynamically the supernova-driven winds in the 
latter. The reasonable agreement (given the large scatter) between 
WF2 and WF2LR suggests that this result is not unduly influenced 
by numerical resolution. 

Although the average galaxy formation efficiency depends 
strongly on feedback, its dependence on halo mass is weak: 77gai 
varies by less than a factor of 2 over the factor of ~ 30 range in 
virial mass spanned by the simulations. In the absence of feedback 
77gai peaks at low masses: feedback is clearly needed to counter the 
high efficiency of gas cooling in low mass haloes. We highlight 
however that the dependence of 77gai on halo mass must become 
significantly stronger for halo masses below those studied here (i.e. 
Afvir < 10"h^^ Mq) in order to successfully reproduce the mea- 
sured faint end of the luminosity function. 



same angular momentum as galaxies in the AGN run. This is strik- 
ing, since their baryonic masses differ on average by a factor of 
~ 10. Feedback affects the mass of a galaxy much more severely 
than its spin: 9 out of 10 baryons in NoF galaxies are missing from 
AGN galaxies, but their specific angular momenta are, on average, 
the same. 

The thin line in the right panel of Fig. |7] shows the results for 
WF2LR. Despite the large scatter, numerical resolution effects are 
clearly noticeable below '--^ 4 x IQ^^h^^ AIq. This corresponds to 
~ 10* particles per halo for WF2LR; extrapolating this to WF2, 
it would mean that our results there should be credible down to 
~ 5 X W^'^h~^ Mq. It would therefore appear as if the main trends 
shown in Fig. |7] are safe from resolution-induced numerical arti- 
facts. 

The angular momentum efficiency peaks in moderate feed- 
back runs (such as WF2) at roughly 50% and its dependence on 
feedback strength is non-monotonic. Despite this apparent com- 
plexity, galaxy ma sses and angular m omenta are actually well cor- 
related. Following ISales etaL j2009l) . we define the mass and an- 
gular momenlum fractions nid and jd, as 



3.7 Galaxy Angular Momentum 



md = r?gal /bar — Mgal/A/vir 



(2) 



The size and rotation speed of galaxy disks place powerful ob- 
servational constraints on galaxy formation models, and are di- 
rectly linked to the angular momentum acquired and retained by the 
baryons that make up the galaxy. We explore this in Fig [7] where 
we show, in the left panel, the specific angular momentum of the 
various galaxy components as a function of halo virial mass. 

As in Fig. [6] dots correspond to the baryonic component in- 
side Tgai for individual systems in run WF2. Although the scatter 
is large (an rms of ~ 0.27 dex), the solid curve, which tracks the 
median j as a function of i\fvir, shows that the specific angular 
momentum scales roughly like j oc M^^'"^. This is the same scal- 
ing found for the dark matter component within rvir (dashed black 
line), and is indeed the expected scaling if the dimensionless halo 
spin parameter, A = J\E\^^^ /GM^^^ , is constant. 

When all baryons within the virial radius are considered, their 
specific angular momentum agrees well with that of the dark mat- 
ter (magenta dotted line). On the other hand, the specific angular 
momentum of central galaxies is, on average, about 50% that of its 
surrounding halo. This fraction, which we refer to as the "angular 
momentum efficiency", rjj — jgai/jvir, appears, on average, to be 
roughly independent of halo mass for WF2 galaxies. 

Interestingly, the gaseous and stellar components of WF2 
galaxies have distinctly different angular momenta. The gas has 2 
to 3 times larger specific angular momentum than the stars, im- 
plying that the radial extent of gaseous disks in these galaxies is 
substantially la rger than that of th e stellar component. This was al- 
ready noted bv lSales et alj ^20091) as a possible way to explain the 
large sizes of the z = 2 star-forming (gaseous ) disks analyzed by 
the SINS survey jForster Schreiber et al.ll20()^ . We shall return to 
this issue in Sec.|4] 



3. 7. 1 Feedback dependence 

The feedback dependence of the angular momentum efficiency, rjj , 
is shown in the right panel of Fig. [7] as a function of virial mass. 
Like the galaxy formation efficiency, ?7gai, the halo mass depen- 
dence of rjj is weak. Unlike T^gai, however, rjj depends only weakly 
on feedback. In a given halo, NoF galaxies have approximately the 



and 



Jd 



(Tvir 



Mgal jgal 
7\</vir jvir 



= '7gal77j/ba 



(3) 



These parameters were introduced by IMo et al] ( Il998h . and have 
become standard fare in semianalytic models of disk galaxy forma- 
tion. 

ISales etaP (120091 noted that jd and rrid correlate well, but in a 
manner different f rom the typical j d = rrid assumption of semiana- 
lytic models (e.g.. lCole et al.l200g) and, perhaps more importantly, 
insensitive to feedback. These authors showed that the simple ex- 
pression 



jd = 9.71 ml (1 - exp[-l/(9.71 rud)]) 



(4) 



provides a good approximation to the results of four OWLS runs 
with supernova-driven winds: WFI, WF2, WF4 and WF2Dec. 

We revisit this result in Fig. [8] where we show the jd-rrid cor- 
relation for all the OWLS runs considered in this paper. The dots 
show individual WF2 galaxies, and are meant to illustrate the typ- 
ical scatter in the relation; the curves trace the median trend of jd 
with rud for the different runs while the black dotted curve outline 
the relation in Eq. |4] Although the AGN and NoF galaxies deviate 
somewhat from the trend outlined in eq.|4] (indicated by the dotted 
thick line), the departures are relatively small and the agreement be- 
tween runs seems remarkable given the extreme range in feedback 
models explored here. 

The bottom panel in Fig. [8] shows the distribution of md for 
four different feedback implementations. Clearly, feedback, at least 
as implemented in our models, affects mostly the baryonic mass of 
galaxies, but largely preserves the link between the spins of haloes 
and galaxies. This link imprints correlations between galaxy mass, 
size, and rotation speed that may be contrasted with observations. 
We turn to this issue next. 
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Figure 7. Left: Specific angular momentum, j, as a function of virial mass. Tlie black dashed line tracks the median j of the dark matter component as a 
function of A/vir. This follows closely the j oc M^^^ correlation expected for systems with constant spin parameter, A. The other symbols, colors, and line 
types are the same as in Fig.|6] Note that the specific angular momentum of all baryons within r^ii is quite similar to that of the halo as a whole (top dotted 
curve). The specific angular momentum of the central galaxy is typically lower than that of the halo; although it correlates well with Afvir, the scatter is large. 
Right: Angular momentum "efficiency", rjj = igai/jvir. as a function of mass for various runs. For clarity, only the median of galaxies in runs NoF, WF2, 
WF2Dec, and AGN, are shown as a function of mass. Note that, unlike r^g^i, the angular momentum efficiency, rij, is a weak function of both mass and of 
feedback. See text for further discussion. 



4 OBSERVATIONAL DIAGNOSTICS 



4.1 Size and Stellar Mass of « = 2 Galaxies 



The feedback-driven trends of galaxy mass and angular momentum 
efficiencies discussed above imprint different relations between tiie 
stellar mass, M^tr, and the size of a galaxy. This is shown in Fig.[9l 
where the various panels compare (for runs NoF, WF2, WF2Dec, 
and AGN) the half-mass radius of the galaxy versus A/str. The pan- 
els on the left show the the half-mass radius of the gas component 
whereas those on the right correspond to the stars. The thick solid 
curve in each panel traces the median trend as a function of Mstr- 
As noted above, simulated galaxies are substantially more extended 
in gas than in stars. 

The simulated galaxies are contrasted with data for the 
large star-forming gas dis ks studied by the SINS survey 
jForster Schreiber et alj|200sl . her eafter FS09), as well as w ith the 
quiescent compact red galaxies of Ivan Dokkum et al. I j2008l. here- 
after vD08). These two datasets probably bracket the extremes in 
the size distribution of massive galaxies at 2: = 2, from the most 
extended to the most compact. 

When feedback is inefficient (e.g., the NoF run) most stars 
form in dense, early-collapsing progenitors that merge later 
on to form the final galaxies. During such mergers the bary- 
onic component transfers angular momentum to the surrounding 
halo, lead ing to the formation of very compa ct massive galaxies 
(iNavarro & Benz 1991; Navarro et al. 1995; N avarro & Steirmiet3 
1 1997 ). The galaxies that result are therefore nearly as compact as 
the quiescent vD08 spheroids, although we note that many of those 
simulated galaxies have half-mass radii even smaller than the grav- 
itational softening of our simulations, so their true sizes are actually 



uncertain. The gaseous component in these simulations is also quite 
compact, with radii rarely matching those of SINS disks. 

Intermediate strength feedback (e.g., the WF2 run) has little 
effect on the most massive galaxies, which are generally as compact 
as the vD08 spheroids. On the other hand, feedback affects more 
strongly less massive systems, leading to a correlation between the 
mass and size of the stellar component where, at the massive end, 
size decreases with increasing mass. This trend runs counter the 
well-established galaxy scaling laws at z=0 (brighter galaxies tend 
to be bigger). The trend is reversed at lower masses and results, 
overall, in systems whose gaseous disks overlap in properties with 
those of galaxies in the SINS survey. 

Increasing the effects of feedback (as in the WF2Dec and 
AGN runs) continues this trend, gradually reducing the mass of 
galaxies and increasing their size at given A/str- This is because 
the more efficient the feedback the more massive the halo inhab- 
ited by a galaxy of given stellar mass. More massive haloes are 
larger and have higher specific angular momenta. Since, as we saw 
above, galaxies generally inherit the specific angular momenta of 
their surrounding haloes, it is possible to have fairly large galaxies 
of modest stellar mass because they actually inhabit large, massive 
haloes. Indeed, many gaseous disks in the AGN run are even more 
extended than the rather extreme examples surveyed by SINS. 



4.2 The Tully-Fisher Relation at z = 2 

The structural diversity of 2: = 2 galaxies discussed above should 
also be manifest in their kinematics. We explore this in Fig. [TO] 
where we plot, as a function of stellar mass, the circular velocity 
estimated for SINS galaxies and for the compact vD08 galaxies. 
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Figure 9. Left: half-mass radius of the gas as a function of stellar mass. Solid black dots in each panel show the results for four of our simulations NoF, WF2, 
WF2Dec and AGN. The thick solid line tracks the median as a function of mass. Open symbols with error bars correspond to the extended star-forming disk 
galaxies from the SINS survey (Forster Schreiber et al. 2009), while the red shaded ellipsoid indicates the area of the plot occupied by the sizes of the compact 
quiescent galaxies from, van D okku m et al. ( 2008). Right: same as be fore, but for t he half-mass radii of the stars. In this case, open symbols and error bars are 
used to indicate the sizes of the compact stellar spheroids from v an Dokkum et a i]f2008), while the shaded blue a rea indicates the region of the plot occupied 
by e xtended gaseous disks fro m the SINS sample. The extended disks reported by SINS jForster Schreiber et al.l 2009) and the compact spheroidal galaxies 
from lvan Dokkum et al ] j2008l) probably bracket the size distribution of massive galaxies at z = 2. The black shaded area indicates the gravitational softening 
of the simulations. Half-mass radii for the gaseous components of simulated galaxies are typically larger than for the stars. Note that the size-stellar mass 
correlation is heavily dependent on feedback. When feedback is very efficient (e.g., WF2Dec/AGN) the size of the gaseous disks increase with stellar mass, a 
correlation that is reversed when feedback efficiency is low. 



For SINS, we use the "maximum" gas rotation speed, as quoted 
by FS09, whereas for vD08 we estimate the circular velocity at the 
effective radius based only on the contribution of the stellar com- 
ponent; i.e., — G'(Afstr/2)/i?cfi. This is clearly a lower limit 
to the circular velocity at that radius, since it neglects the possible 
contributions of dark matter and gas components. We note this in 
Fig.[TO]by small arrows on the vD08 data points (open triangles). 

It is clear from this rendition of the data that the two popula- 
tions of z — 2 galaxies follow very different Tully-Fisher relations. 
At given stellar mass, the compact galaxies are expected to have 
circular velocities at least twice higher than SI NS disks. Although 
kinem atic data for such galaxies is scarce, Ivan Dokkum et all 
( l2009h report a preliminary measurement of the velocity disper- 
sion of one of these galaxies. The high velocity dispersion reported, 
~ 510 km/s, agrees with this interpretation. 

The circular velocity of the simulated galaxies is measured at 
the half-mass radius of the stellar (red solid curve) or the gaseous 
(blue dashed curve) component, respectively. The comparison be- 
tween simulations and observations yields similar conclusions as in 
the previous subsection. 

Inefficient or absent feedback (e.g., NoF) yields galaxies that 
are more concentrated than the SINS disks, and therefore have, 
at given stellar mass, typically higher circular velocities. Forming 
large, extended disks is difficult in the absence of efficient feed- 
back. By contrast, accounting for the compact spheroids studied by 
vD08 is relatively easy. 

In the case of AGN or WF2Dec, the most efficient feedback 
schemes explored in Fig. [TO] many simulated galaxies are as spa- 



tially extended as the SINS disks, and the good agreement extends 
to the Tully-Fisher relation for those galaxies. Few very massive 
galaxies form as a result of the efficient feedback, and very few of 
those that form are as compact as those in the vDOS sample. 

More moderate feedback choices give intermediate results. 
We consider it encouraging that some galaxies in the WF2 runs 
overlap with both SINS and vD08 in Figs.|9]and [TO] If these mod- 
els are correct, then there should be a sizable population of galaxies 
at z — 2 with properties intermediate to the SINS disks and vD08. 

To summarize, the results shown in Figs. |9] and 1 101 indicate 
that neither the extreme compact sizes of massive spheroids nor 
the large spatial extent of star-forming disks at z — 2 pose insur- 
mountable challenges to the standard paradigm. Indeed, it is possi- 
ble, with adjustments to the feedback algorithm, to reproduce either 
population without resorting to unusual halo spin or halo formation 
histories. 

At the same time, reproducing the striking diversity in the ob- 
served sizes and masses of z = 2 galaxies with a single feedback 
recipe might be challenging, but we are encouraged by the large 
scatter in the properties of simulated galaxies at given stellar mass 
that arises naturally in any feedback model. The relative abundance 
of either population is still poorly constrained observationally, and 
our small simulation box might not be adequate to study or search 
for rare, extreme populations. Improved observational constraints 
on the relative abundance of extended vs compact galaxies and a 
better characterization of the "average" population of 2: = 2 galax- 
ies will certainly help to constrain which feedback implementation 
gives results that agree best with observation. 
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Figure 8. Top: The angulai' momentum fraction = Jgai/>^vir the 
galaxy mass fraction, = Mgai/Afvir- I3ots con'espond to individual 
galaxies in WF2 and ai'e meant to illustrate the scatter; the median trend 
is traced by the black solid line. Other curves are analogous, but for each 
fe edback mod e l anal yzed here. The black dotted curve is the fit proposed 
bv lSales et alj j2009t) . The straight hne labeled = corresponds to 
the commonly-adopted assumption that the specific angular momentum of 
a galaxy equals that of its surrounding halo. Bottom: Distribution of galaxy 
mass fraction, m^, for four different runs spanning the range of feedback 
strengths of our simulations: NoF, WF2, WF2Dec, and AGN. 



Figure 10. The stellar mass-circular velocity (Tully-Fisher) relation for 
z = 2 galaxies identified in runs with four different feedback implemen- 
tations. Symbols are as in Fig. |9] The median circular velocity measured 
at the stellar half-mass radius is shown by the solid red line. Vertical lines 
show the 25-75 percentiles of the distribution. We also show this relation 
when the circular velocity is measured at the half-mass radius of the star 
forming gas (dashed blue curve). Open circles and triangles show the obser- 
vational determinations for di sks a nd compact galaxie s at z= 2 taken from 
iForster Schreiber et alj j2009h and Ivan Dokkum et al.l <2008l) . For the lat- 
ter we assign velocities by neglecting the dark matter distribution; i.e., we 
assume = G(Mstr/2)/i?cff. which constitutes a lower limit to the 
true circular velocity. This is indicated by the horizontal arrows in each 
panel. The thick dotted line is the Bell & de Jond l200ll) relation for late- 
type galaxies at z=0 corrected to a Chabrier IMF. 



4.3 Disks and Mergers at z = 2 

Another interesting constraint is provided by the frequency of sys- 
tems actively forming stars in rotati onally-supported disks. Before 
surveys such as SINS and OSIRIS jporster Schreiber et alj[2009l : 
iLaw et alj2009l and references therein) started to resolve the kine- 
matics of star-forming galaxies at high z, it had been common- 
place to assume that systems where star formation was progress- 
ing in earnest would almost invariably be ongoing major merg- 
ers. It is now clear, however, that at least about one third of the 
galaxies surveyed by SINS and OSIRIS are forming stars in rela- 
tively quiescent disks rather than ongoing mergers with disturbed 
and transient k inematics (for an alternative view, however, see 
[Robertson & Bullock,2008i) . 

We use Krot , the simple measure of the importance of ordered 
rotation introduced in Sec. 13.41 to explore this issue in our simu- 
lations. When most of the gas is in a rotationally-supported disk, 
the parameter Urot should approach unity. Fig.[TT]enables a visual 
calibration of this parameter by showing edge-on projections of 12 
galaxies arranged by the value of Krot of the central galaxy (in this 
case only the star-forming gas is used to compute Krot)- Figure [TT] 
shows an image of the projected gas density within a sphere of ra- 
dius l.Srgai. Thin, extended disks are the norm when Krot<;0.75. 
Ongoing mergers typically have Krot;$0.5; those with intermediate 
values of Krot have disturbed morphologies, and tend to be late- 
stage mergers or systems where accretion is ongoing but minor 



Using this simple measure, the fraction of ongoing mergers 
vs quiescent disks may be readily estimated, and is shown in the 
top panel of Fig. [12] for the case of WF2. The distribution of Krot 
for all WF2 galaxies is shown by the top histogram; the shaded 
histogram is for the same run, but reducing the sample of galax- 
ies to one half by selecting only those in haloes more massive than 
2x10^^ h"^ Mq. Encouragingly, the shape of the two histograms is 
quite similar This is further confirmed by the distribution of Krot in 
WF2LR galaxies (for Mvir > 2 x lO^/i^H/e) which is shown as 
the thin solid line in the bottom panel of Fig. [12]. The good agree- 
ment between WF2 and WF2LR indicates that numerical resolution 
effects are unlikely to compromise our conclusions. 

According to the definition above, about 45% of WF2 galaxies 
are reasonably quiescent star-forming disks, and only about 20% 
are ongoing major mergers. These fractions are similar for WF2 
and WFl, and seem consistent with the observational data quoted 
above. 

For the run without feedback, NoF, over ~ 75% of the galax- 
ies are classified as disks. This is because, in the absence of feed- 
back, the gas cools and flows unimpeded to the center, where it 
settles into disks and forms stars profusely. These disks are, how- 
ever, quite small (see Fig.|9]l. The absence of effective feedback al- 
lows the gas to remain undisturbed in such disks, which are quickly 
reconstituted after me rgers (see, e.g.. lSDringel & Hernauisl|[2005l ; 
[Robertson et al.[|2006l) . At the other extreme, only 5% of all galax- 
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ies in the AGN run, and ^ 20% of those in the WF2Dec run, would 
be classified as disks according to this criterion. 

Strong feedback-driven winds can clearly disturb quiescent 
disk morphologies, and their kinematic effects may be difficult to 
disentangle from those of ongoing mergers. It remains to be seen 
whether a simple feedback model can account for both the ob- 
served frequency of galaxies with disk-like kinematics as well as 
the mounting evide nce for large scale galactic outflows z ^ 2 
dSteidel et al.ll2010h . 



5 SUMMARY AND CONCLUSIONS 

We study the effects of various feedback implementations on the 
structure and morphology of simulated galaxies at 2 = 2. Our 
analysis uses nine runs from the Overwhelmingly Large Simula- 
tions (OWLS) project, and probe a variety of possible feedback im- 
plementations, from "no feedback" to supernova-driven wind feed- 
back to strong outflows aided by the contribution from AGNs. Ex- 
cept for the no-feedback and AGN-feedback cases, all other runs 
assume that the same amount of feedback energy (per mass of stars 
formed) is devolved by supernovae to the interstellar medium: the 
main difference is how this energy is coupled to the medium, which 
in turn determines the overall effectiveness of the feedback. 

Each run follows the evolution of the same 25 Mpc box 
up to = 2, with 512^ dark matter particles and 512"^ parti- 
cles for the baryonic component. All other simulation parame- 
ters (star formation algorithm, stellar initial mass function, etc) 
are kept constant, so any differences between runs may be traced 
solely to feedback. In total, we analyze for each run ~ 150 galax- 
ies formed at the centers of haloes with virial mass in the range 
10" Mq < A/vir < 3 X lO^^/i"^ Mq. Our main results may 
be summarized as follows. 

• Varying the feedback implementation can lead to dramatic dif- 
ferences in the mass of galaxies formed in a given dark matter halo. 
The galaxy formation efficiency, 77gai — Afgai/(/barA/vir), varies 
by roughly an order of magnitude when comparing the no-feedback 
run (NoF, where r^gai ~ 0.5) to the AGN-l-supernova feedback run 
(AGN, where r/gai ~ 0.05), the two extremes probed by our simu- 
lations. 

• The ability of feedback to regulate the efficiency of galaxy 
formation in haloes of different mass varies according to the details 
of the adopted numerical implementation of the feedback. Weak 
or ineffective feedback leads to a decrease in galaxy formation ef- 
ficiency with mass, whereas strong feedback curtails preferentially 
the formation of galaxies in low-mass haloes. The mass dependence 
is, however, modest, with variations in r^gai of less than a factor of 
~ 2 over the (factor of ^ 30) mass range spanned by haloes in our 
sample. 

• Feedback results in strong correlations between galaxy mass 
and angular momentum. This leaves an imprint on galaxy mor- 
phologies and on the scaling laws relating mass, size, and circular 
velocity. 

• Weak feedback minimizes disturbances to the settling of gas 
in rotationally-supported structures, and favors the formation and 
survival of quiescent gaseous disks. However, weak feedback also 
allows much of the gas to form stars early in dense protogalactic 
clumps that are later disrupted in mergers as the final galaxy as- 
sembles. Such mergers also transfer angular momentum from the 
baryons to the halo. The net result is a predominance of dense. 



spheroid-dominated stellar components and a scarcity of spatially- 
extended star-forming disks. 

• Strong feedback, on the other hand, promotes the formation of 
large, extended galaxies. Indeed, the more efficient the feedback the 
more massive (and therefore, larger) the halo inhabited by a galaxy 
of given stellar mass. It is thus possible to have fairly large galaxies 
of modest stellar mass because, when feedback is strong, they in- 
habit large, massive haloes. The size, mass, and rotation speeds of 
these extended galaxies compare favorably with those reported by 
the SINS survey. This, however, comes at the expense of inhibiting 
the survival of rotationally-supported disks of quiescent kinematics 
and of preventing the formation of compact stellar spheroids. 

• Moderate-feedback runs result in galaxies that follow scaling 
laws that are intermediate between large s tar-forming disks, such 
as tho se studied by the SINS collaboration JForster Schreiber et al.l 
'200P), and the compact, quiescent early-type systems analyzed by 
van Dokkum et al. (2008). Disk-like morphologies in both gas and 
stars are common in these runs, in numbers that appear commensu- 
rate with current constraints. 

Although far from definitive, the results outlined above are en- 
couraging. Properly calibrated, simple feedback recipes such as the 
ones we explore here seem able to produce galaxies with properties 
in broad agreement with observation. One should be aware, how- 
ever, of the numerical sensitivity of the results to details of feedback 
implementation. Nevertheless, if developed in step with observa- 
tional progress in the characterization of the high-redshift galaxy 
population, simulations are likely to become more and more reli- 
able tools, useful when trying to make sense of the striking diver- 
sity of high-z galaxies in terms of the current paradigm of structure 
formation. 
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